ARL-TN.58 


AR-008^ 


DEPARTMENT  OF  DEFENCE 

DEFENCE  SCIENCE  AND  TECHNOLOGY  ORGANISATION 
AERONAUTICAL  RESEARCH  LABORATORY 

MELBOURNE,  VICTORIA 


Technical  Note  58 

STRESS  ANALYSIS  OF  A  SECONDARY-BENDING  SPECIMEN 


by 

RX.  EVANS 
M.  HELLER 


1 


Approved  for  public  release 


©  COMMONWEALTH  OF  AUSTRALIA  1993 

NOVEMBER  1993 


W^T*~24*“0T4 


^  94-02036 


This  work  is  copyright  Apsit  from  any  uss  as  pormitted  under  ths 
Copyright  Act  1968,  no  port  may  bo  rsproducad  by  any  process  vrithout 
prior  written  permission  from  the  Austraiian  Government  Pubiishing 
Services.  Requests  end  enquiries  concerning  reproduction  and  rights 
shouid  be  addressed  to  the  Manager,  Commonwealth  information 
Services,  Austraiian  Government  Pubiishing  Services,  GPO  Box  84, 
Canberra  ACT  2601. 


THE  UNTTEO  STATES  NATIONAL 
TECHNICAL  INFORMATION  SERVICE 
IS  AUTHORREO  TO 

rerroouce  and  sat  this  report 


DKClAQIll  NOTICE 


TfflS  DOCUMENT  IS  BEST 
QUALITY  AVAILABLE.  THE  COPY 
FURNISHED  TO  DTIC  CONTAINED 
A  SIGNIFICANT  NUMBER  OF 
COLOR  PAGES  WHICH  DO  NOT 
REPRODUCE  LEGIBLY  ON  BLACK 
AND  WHITE  MICROFICHE. 


AR-OOg-404 


DEPARTMENT  OF  DEFENCE 

DEFENCE  SCIENCE  AND  TECHNOLOGY  ORGANISATION 
AERONAUTICAL  RESEARCH  LABORATORY 


Technical  Note  58 


STRESS  ANALYSIS  OF  A  SECONDARY-BENDING  SPECIMEN 


RX.  EVANS 
M.  HELLER 

SUMMARY 


This  Note  describes  a  two-dimensional  finite-element  elastic  analysis  of  a  uniaxialfy-loaded 
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L  INTRODUCTION 


A  test  {sogram  has  recently  been  undertaken  at  ARL  to  invesdgaie  the  effect  of  secondary 
bending  on  die  fatigue  life  of  jc^ts  with  and  without  hole  cold-working  and  inteifeienoe-fit 
bolts  [1].  SectMidaiy  bending  is  caused  by  an  eccentricity  of  the  applied  load  relative  to  die 
neutral  axis  of  a  structural  ccunponent,  and  results  in  the  superposition  of  bending  ^sses 
on  axial  stresses.  Hiis  Note  describes  a  preliniinary  two-dimensional  finite-element  elastic 
analysis  oi  a  uniaxially-loaded  bolted  seomdary-bending  qiecimen  which  was  c(»ducted  to 
infocmatkm  rdevant  to  the  ARL  fatigue  testing  progranL  Thermoelastic  tests  were 
amducted,  and  the  experimental  results  are  compared  widi  the  finite-element-model  analysis. 
Experimental  information  from  the  United  Kingdom  on  similar  qiedmens  is  also  examined. 


2.  SPECIMEN  DETAILS 

The  secondary-bending  specimen  analysed  in  this  paper  is  depicted  in  Hgure  1,  and  is  a 
slighdy  modified  version  of  the  U.K.  Q-joint  [2].  It  is  composed  of  three  main  members 
(two  S  mm  diick  components,  and  one  3  mm  thick  component)  fastened  together  with  four 
bdts  of  6.30  nan  diameter.  A  2  mm  and  a  5  mm  spacer  conqilete  the  specimen.  (The 
modification  with  respect  to  the  UJC.  (^jdnt  is  an  extension  of  the  gripping  pmtion  of  the 
qiecimen  from  60  mm  to  7S  mm  to  allow  easier  gripping  in  ARL  test  mattes).  The 
material  is  a  high  strength/weight  aluminium  alloy,  namely  707S-T6S1,  commonly  used  in 
the  F/A-18  fighter  aircraft,  with  die  elastic  properties  E  =  71  GPa  and  v  =  0.33. 


3.  NUMERICAL  ANALYSIS 

A  total  eight  models  of  differing  degrees  of  rigidity  were  analysed.  The  three  discussed  in 
this  report  cover  the  range  that  were  analysed  and  vary  from  a  totally  rigid  joint  to  a  mudi 
more  flexiUe  and  realistic  joint  The  models  are  two-dimensional,  assuming  plane  stress, 
while  the  actual  specimen  is  of  course  three-dimensitxial,  pertinently  in  the  region  of  the 
bolts.  On  comparing  Hgure  1  of  the  specimen,  to  Hgure  2  of  the  finite-element-noodel 
getnnetry,  it  can  be  seen  ^t:  (i)  the  actual  gripping  portion  of  the  specimen  is  not  rrxxlelled 
(however,  the  restraints  are  correcdy  represented),  and  (ii)  actual  bolts/holes  are  not 
included  in  the  noodeL  Model  1  incorporates  rigid  joints  between  the  components.  Model  2 
and  Model  3  use  'repeated  freedoms'  ami  'gap  elements',  at  the  interface  betweoi 
components  in  an  effort  to  model  the  bolts/holes,  and  hence,  the  bending  behaviour  of  the 
bolted  joint  Model  3  is  andlar  to  Model  2,  the  difference  being  the  number  of  repeated- 
fieedmirs  pairs  and  hence  gap  elements. 

Repeated  fipeedocos  are  employed  to  restrict  linked  nodes  (pairs)  to  move  identically,  hence, 
they  always  remain  a  fixed  distance  apart  (0.001  mm,  for  Models  2  and  3).  They  noodel  the 
areas  of  the  specimois  where  the  components  are,  effectively,  cotmected  by  the  ccurqnessive 
stresses  produced  by  the  toiqued  bolts.  Gap  elements  are  also  used  to  restrict  the  movement 
of  linked  nodes,  however,  in  a  limited  way.  They  enable  the  adjacent  components  to 
separate  and  slide  but  not  to  ovedap,  and  thus,  model  the  areas  of  the  specimen  adjacent  to 
the  bdts,  which  are  out  of  the  influence  of  the  compressive  stresses  of  tte  bolts. 

The  increase  in  width  of  the  ends  of  the  specimen  is  accounted  for  in  the  two-dimenskxud 
noodels  by  a  change  in  depth  of  the  relevant  elements.  The  increase  in  width  of  the 
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exmsioQs  is  modelled  by  diree  rectangles;  the  total  area  of  the  extensions  oi  the  model 
being  equivalent  to  that  of  die  extensions  of  die  specimen.  The  elements  of  the  controlling 
and  test  sections  have  a  depth  of  51  mm.  The  depth  of  the  elements  at  the  thinnest  end  at 
the  extension  is  S1.5  mm  for  a  length  of  26  mm,  the  depth  then  increases  to  S9.S  mm  for 
13  mm,  and  lasdy,  to  72  mm  for  13  mm  at  the  end  of  the  extension;  the  total  length  of  each 
extension  is  equal  to  the  sum  of  26, 13  and  13  (52  mm).  The  local  neutral  axis  (ynlirection) 
positions  relative  to  that  for  the  tqiplied  load  are:  -1.0  mm  for  the  left-hand  extensitni, 
->-U  mm  for  die  test  section,  0  mm  for  the  controlling  section,  and  40.625  mm  for  the  right- 
hand  extension. 

In  all  diree  models  a  uniform  tensile  stress  is  tqiplied  ro  the  right-hand  end  of  the  model.  The 
left-hand  aid  is  restrained  from  moving  horizontally  and  vertically,  and  all  end  nodes  are 
restricted  to  moving  horizontally  together  (Figure  2(a)).  These  restraints  model  the 
specimen  as  if  it  is  gripped  in  a  fatigue  test  machine.  The  resultant  finite  element  mesh 
(Figure  2(b))  ccxisists  of  582  eight-noded  isoparametric  quadrilateral  elements.  Model  1  has 
1957  nodes  and  Models  2  and  3  have  2035  nodes.  Listings  of  the  finite  element  data  files 
for  Model  1  and  Model  3  are  given  in  Appendices  A  and  B,  respectively.  To  enable 
cooqiarison  between  the  theoretical  and  experimental  stress  distribution  data,  the  bulk-stress 
(Ox  4-  Oy)  contour  plots  for  each  of  the  models,  under  an  applied  tensile  stress  of  16.4  MPa 
are  given  in  Figures  3(a)-(c).  The  stress  value  of  16.4  MPa  is  equivalent  to  the  tnaxiinmn 
tensile  stress  tilled  to  the  experimental  specimen.  Hgure  4  is  the  displacement  diagram 
obtained  from  Model  3. 

The  PAFEC  (Program  for  Automatic  Hnite  Element  Calculations)  suite  of  programs  was 
used  on  an  'Apollo  DN1(X)00*  machine  with  the  'Apollo  Domain/OS  Version  SR10.3' 
operating  system  for  this  work. 

3.1  Model  1 

Model  1  has  'welded'  joints.  The  joints  between  the  three  components  are  rigid,  i.e.,  the 
specimen  is  modelled  as  a  single  machined  article. 

3.2  Model  2 

Model  2  uses  repeated  freedoms  and  gap  elements  to  join  the  three  separate  components.  A 
gtq>  of  0.001  mm  between  the  components  is  introduced.  The  repeated  freedoms  join  pairs 
of  nodes  on  adjacent  components  at  the  bolt  sites  and  extend  either  side  of  the  'bolt' 
diameters  to  the  next  and  following  comer  nodes  (i.e.,  a  further  two  elements  on  each  side 
frtxn  the  edges  of  the  bolts).  A  total  of  forty-one  (41)  separate  pairs  of  nodes  are  forced  to 
move  as  angle  nodes  by  using  repeated  freedoms,  joining  the  components  togedier  in  three 
locations.  This  simulates  the  Yorcing  together'  of  surfaces  of  the  bolted  joint  of  the 
^lecimen.  Gap  elements  are  used  to  link  the  components  at  all  nodes  either  side  of  die  'bolts' 
(thirty-seven  (37)  separate  pairs). 

33  Model  3 

Model  2  was  found  to  be  too  rigid  in  the  vicinity  of  the  bolts,  thus  Model  3  employs  less 
repeated  freedoms  and  more  gap  elements  than  Model  2.  The  repeated  freedoms  in  Model  3 
are  used  at  the  bolt  sites  only.  A  total  of  seventeen  (17)  separate  pairs  of  nodes  are  forced 
to  move  as  angle  nodes  using  repeated  freedoms,  and  sixty-one  (61)  separate  pairs  of  nodes 
are  linked  by  gap  elements. 
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4.  EXPERIMENTAL  ANALYSIS 

SPATE  (Stress  Pattern  Analysis  measurement  of  Thermal  Emission)  testing  was 
perfinmed  on  a  seocmdary-bending  ^)ecimen  had  cold-expanded  boh  holes  and 
interference-fit  bolts  (qwdmen  SQ7S).  The  SPATE  8(XX)  system  used  in  diis  experimemal 
analysis,  detects  emitt^  diennal  energy  fixnn  a  q)ecimen  under  cyclic  stress.  It  scans  a 
preselected  area  and  an  uncalibrated  stress  contour  plot  of  (Gx  Gy)  is  obtained  point  by 
point 

A  ctMistant  ahemating  load  of  0.S  -  10.5  kN  (0.78  -  16.4  MPa)  was  tilled  and  the 
following  scans  taken,  with  'front',  iMck',  'right'  and  'upper*  orientations  as  per  Figure  1 : 

1  Right  side  scan,  frt»n  front  angle  of  15^  off  axis,  with  SPATE  detector-head  56  cm 
from  ^tecimen  (resolution  5), 

2  Ri^t  side  scan,  from  rearward  angle  of  15**  off  axis,  widi  SPATE  detecttx'-head  66  cm 
from  specimen  (resolution  5), 

3  Frcmt  scan  of  upper  bolts,  with  SPATE  detector-head  31  cm  from  specinoen 
(resolution  4),  and 

4  Back  scan  of  upper  bolts,  with  SPATE  detector-head  41  cm  from  specimen 
(resolution  4). 

Scans  1  and  2  can  be  compared  with  finite  element  results  once  the  correqxMiding  area  is 
outlined  (Hgures  5,  6  and  3).  Scans  3  and  4  (Figures  7  and  8)  highlight  the  three- 
dimensional  nature  of  the  stress  distribution  in  the  specimen. 

5.  DISCUSSION  OF  RESULTS 

A  sketch  of  the  joint  sectitxi  of  the  specimen  is  incorporated  for  clarity  in  Hgure  5  (scan  1) 
and  Figure  6  (scan  2).  Hgure  3(d)  is  a  side  view  only  of  Figure  5  which  enables  comparison 
between  the  numerical  stress  contour  plots  of  Figures  3(a)-(c)  and  the  experimental  data. 
Information  regarding  the  secondary-bending  ratio  [3]  is  included  in  Appendix  C. 

5.1  Numerical  Analysis 

In  ctxnparing  the  bulk-stress  (Ox  +  Oy)  contour  plots  of  Model  1,  Model  2  and  Model  3  (in 
Hgures  3(aXc)),  it  can  be  seen  in  all  three  that  the: 

(i)  stress  flows  through  each  component, 

(ii)  main  component  carries  the  highest  average  stress, 

Qii)  stress  decreases  as  y  decreases  in  the  extension  of  the  main  component, 

(iv)  stress  decreases  as  y  decreases  and  then  changes  such  that  stress  increases  as 
y  decreases  in  the  test  section  of  the  main  component, 

(v)  stress  increases  as  y  decreases  in  the  extensions  of  the  secondary  components,  and 

(vi)  'free  ends'  of  the  ctHnponents  carry  a  very  low  (or  negative)stress. 

Hgure  3(a)  also  highlights  the  sharp-edge  stress  raisers.  Hgures  3(b)  and  3(c)  show  the 
singularities  caused  by  a  change  in  stiffness  at  the  interf^  of  Ae  repeated  freedoms  arxl  gq> 
elements.  They  also  indicate  that  the  5  mm  diick  secondary  component  carries  the  least 
stress. 

The  displacemoit  diagram  of  Hgure  4  confirms  the  occurrence  of  points  (iii),  (iv)  and  (v) 
above,  i.e.,  bending.  The  noodels  bend  similariy  to  Hgure  4  due  to  the  neutral  axis  of;  the 
left-hand  extensitm  being  lower  by  1.0  mm,  the  test  section  being  higho’  by  1.5  mm,  and  the 


ri^-hand  extrasion  being  higher  by  0.625  mm,  Uian  that  of  the  rested  die  model  Wheadie 
qiecimen  experiences  an  ^iplied  tensile  stress,  these  three  areas,  due  to  end  constraints,  try 
to  align  dwir  neutral  axes  with  that  of  the  gripped  secdons  Oocated  at  y  -6.5  mm). 

5JZ  Experimental  Analysis 

Hgures  5  and  6  are  SPATE  scans  1  and  2  respectively,  with  an  over-sketch  of  the  qiecimen 
fcH*  clarity.  The  presence  of  the  SPATE  signal  outside  the  boundary  of  the  specimen  is  due 
to  limitations  in  the  depth  of  field  of  the  infrared  detector.  These  two  figures  clearly  show: 

Q)  the  stress  distribution  through  components, 

Oi)  that  the  main  conqionent  experiences  the  highest  average  stress, 

0ii)  that  die  main  component  carries  a  high  stress  in  the  front  face  under  the  upper  set  of 
bdts, 

(iv)  that  the  thin  secondary  component  takes  the  next  highest  average  stress, 

(v)  that  the  thick  secondruy  component  experiences  the  lowest  average  stress, 

(vi)  the  stress  gradient  across  the  thickness  of  each  component,  and 

(vii) that  the  free  ends'  of  the  components  carry  minimal  stress. 

Hgures  7  and  8  indicate  the  following,  where  the  'upper'  and  Tower*  set  of  bolts  are  as  shown 
inHgure  1: 

(i)  die  three-dimensional  nature  of  the  stress  field  in  die  vicinity  of  the  bolt  holes, 

Oi)  that  the  main  emnponent  carries  a  maximum  stress  under  dw  upper  set  of  bdts, 

Oii)  diat  the  thin  compcaient  experiences  minimal  stress  above  the  lower  set  of  bolts, 

(iv)  that  the  main  component  carries  a  greater  stress  on  its  back  (compared  to  its  front),  and 

(v)  that  the  thick  secondary  component  experiences  a  small  amount  of  stress  immediately 
under  the  upper  set  of  bolts. 

Therefore,  the  SPATE  scans  provide  infexmation  regarding  the  stress  distribution,  including 
die  bending  characteristics  of  the  specimen. 

S3  Comparison  of  Experimental  and  Numerical  Analyses 

The  general  trend  of  the  specimen's  load-transfer  characteristics  can  be  obtained  from  the 
finite  elonent  models  as  seen  in  Section  5.1.  The  bending  behaviour  of  the  experimental 
specimen  is  also  consistent  with  the  displacement  diagram  of  Model  3  (Figure  4). 

However,  local  values  gained  from  the  numerical  analysis  do  not  correlate  well  with  the 
experimentally  acquired  data  at  certain  regions.  Figure  9  (and  Appendix  C).  The  capability 
of  the  finite  element  models  to  represent  an  actual  specimen  is  summarised  in  Hgure  9. 
Average  values  of  bulk  stress  were  obtained  at  the  6ve  designated  points  (A,  B,  C,  D,  and 
E)  from  each  of  the  three  finite  element  naodels  and  from  the  SPATE  scan  number  1.  Using 
die  axis  system  of  Figure  2(a),  the  coordinates  of  die  five  locations  are  as  fdlows: 

Point  coordinates 
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Points  A  to  E  are  located  in  the  centre  of  the  thickness  of  the  components,  and  the  bulk- 
stress  value  of  each  point  was  determined  by  taking  the  average  of  it  plus  its  four  adjacoit 
data  pc^ts,  where  the  data  points,  for  the  numerical  models,  are  elemoits,  and  fm  the 
erqterimental  analysis,  are  SPATE  measuring  points.  For  each  case  (three  numaical  and  one 
cjqKrimental),  the  value  of  bulk  stress  at  points  B.  C,  D,  and  E,  were  compared  to  the  value 
at  A.  The  relative  bulk  stress  at  points  B  and  C  show  compliance  for  all  cases,  however,  the 
stresses  of  the  numerical  models  differ  from  the  experimental  data  at  point  E,  and  eqrecially 
at  point  D.  Model  1,  the  "welded  joint'  model,  has  relative  bulk-stress  values  at  points  D  and 
E  which  have  much  higher  values  than  those  for  SPATE,  indicating  that  a  much  greater 
relative  stress  is  experienced  by  the  "free  ends'  of  the  components.  Models  2  and  3  show  the 
reverse;  the  free  ends'  sustain  negligible  or  negative  stress  compared  to  SPATE.  Therefore, 
a  duee-dimensional  numerical  model  which  displays  characteristics  between  those  of 
Model  1,  and  Models  2  and  3,  is  desired. 

As  expected,  the  di^ering  finite  element  approaches  bounded  the  true  solution.  It  is  believed 
that  a  three-dimensional  model  is  necessary  to  obtain  good  agreement  between  numerical 
and  experimental  results  at  any  location. 

5.4  Comparison  of  U.K.  and  Numerical  Data 

Experimental  secondary-bending  ratio  data  of  similar  specimens  from  the  U.K.  [2]  were 
examined.  Values  of  secondary-bending  ratio  were  calculated  from  the  finite  element 
models  and  compared  to  the  experimental  infcmnation  from  the  U.K.  specimens  (for  noore 
detail  refer  to  Appendix  C).  llie  secondary-bending  ratio  is  defined  as  'the  ratio  of  the 
bending  strain  and  the  axial  strain  at  the  position  under  consideration'  [3],  and  the  values 
from  the  various  models,  at  a  stress  equal  to  that  used  for  the  SPATE  measurements,  are  as 
follows; 


Secondary-bending  ratio  at  an  applied  stress  of  16.4  MPa 


Model 

No. 

Secondary-bending 

ratio 

1 

-0.3616 

2 

-0.2186 

3 

0.5446 

Model  3  predicts  the  correct  sign  of  the  secondary-bending  ratio,  however,  the  values  of 
secondary-bending  ratio  for  Model  3  are  unduly  increased  by  stress  concentrations  at  the 
nearby  repeated-freedom/gap-element  interface.  The  values  of  secondary-bending  ratio  for 
all  three  numerical  models  remain  virtually  constant  with  change  in  applied  stress,  in  contrast 
to  the  experin^ntal  data,  where  the  values  of  secondary-bending  ratio  increase  with  increase 
in  stress.  The  experimental  U.K.  values  of  secondary-bending  ratio  range  from  0.13  at  an 
t^plied  stress  of  17  MPa,  to  0.494  at  the  maximum  applied  stress  of  104.7  MPa.  As 
menti(xied  above,  it  is  believed  that  a  three-dimensional  model  is  necessary  to  obtain  good 
agreement  between  numerical  and  experimental  results  at  specific  locaticms. 


* 
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€.  CONCLUSIONS 

1.  This  Note  has  shown  that  the  general  trends  of  the  load-transfer  characteristics  of  the 
secondaiy-bending  ^>edinen  can  be  obtained  from  two^iimensional  finite  element  models. 

2.  Stresses  fimn  specific  regions  of  a  two-dimensional  finite  element  model  do  not 
satisfBctorily  con^Moe  to  tire  behaviour  of  an  actual  secondary-bending  bolted  joint 
qredmen.  Difficulty  was  experienced  in  obtaining  the  correct  lo^  transfer  through  the 
l»lts.  Model  1  resulted  in  too  much  load  transfer  and  Models  2  and  3  in  too  little.  A 
dreoretical  model  which  gives  results  between  a  welded  model  and  one  using  repeated 
fieedoms  and  gap  elements  is  desired. 

3.  The  values  of  secondary-bending  ratio  for  the  two-dimensional  models  did  not  agree 
with  those  obtained  experimentally  (from  the  U.K.).  Of  the  three  models.  Model  3  gave  the 
mcret  reasonable  results.  The  experimental  values  of  secondary-bending  ratio  were 
calculated  using  strain-gauge  data.  The  location  of  these  strain  gauges  was  very  close  to  the 
boh  holes.  Qearly,  three-dimensional  effects  dominate  in  the  region  of  bolt  holes,  thus,  a 
two^hmensional  model  is  inadequate.  It  is  believed  that  a  thiee-dinrensional  model  to 
represent  the  secondary-bending  specimen  is  essential  for  the  generation  of  more  consistent 
results. 
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FIGURE  2:  Two-dimensional  finite  element  model 


(c):  Contour  plot  of  bulk  stress  (Gx  +  Cy)  for  Model  3,  with  applied  stress  of  16.4  MPa 


units 


(d):  SPATE  scan  number  1  at  an  applied  stress  of  16.4  MPa 


P’lGURE  3:  Comparison  of  numerical  and  experimental  bulk-stress  distributions 


FIGURE  4:  Displacement  diagram  frcnn  Model  3  under  a  tensile  stress 
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FIGURE  5;  Right  side  SPATE  scan  from  front  angle  (scan  1) 


Right  side  SPATE  scan  from  back  angle  (scan  2) 


FIGURE  7:  SPATE  scan  of  front  of  specimen  (scan  3) 


FIGURE  8:  SPATE  scan  of  back  of  specimen  (scan  4) 
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APPENDIX  A:  PAFEC  data  file  for  Model  1 


THE  SECONDARY  BENDING  SPECIMEN  IS  MODELLED  AS  A  SOLID  COMPONENT 
-  A  'WELDED'  JOINT.  THE  LHS  IS  FIXED,  THE  RHS  IS  RESTRAINED  IN 
THE  VERTICAL  PLANE  ONLY  AND  THE  LOAD  APPLIED  IN  THE  HORIZONTAL 
PLANE,  THE  'WAISTING'  OF  THE  SPECIMEN  IS  MODELLED  IN  THREE 
RECTANGULAR  STEPS.  NOTE:  Y (DATA  FILE)  »  -Y (FIGURES  IN  REPORT) 


C  DATA  FOR  THE  2D  ANALYSIS  OF  A  SECONDARY  BENDING  SPECIMEN 
C 
C 
C 
C 
C 

c 
c 

CONTROL 
FULL. CONTROL 
TOLERANCE=lE-4 
PHASE=1 
PHASE*=2 
PHASE=4 
PHASE=6 
PHASE=7 
PHASE=9 
USE . STREN 
CONTROL.  END 
C 
C 

NODES 

NODE .  NUMBER 
1 


X 

0 

12.5 
25.0 

37.5 
50.0 
102.0 


PAFBLOCKS 

ELEMENT . TyPE=3  6210 
PROPERIES=ll 
BLOCK .  GROUP , 

1 
2 

3 

4 

5 

8 

9 

10 
11 
12 


N1 

1 

2 

3 

4 
1 
2 

3 

4 

3 

4 


1 
2 

3 

4 

7 

8 

9 

10 
11 
12 
C 
C 

PAFBLOCKS 

ELEMENT . TYPE=3  6210 
PROPERTIES=12 
BLOCK.  GROUP.  N1 

14  14  10 

15  15  10 

16  16  10 
C 

C 

PAFBLOCKS 

ELEMENT . TYPE=3  6210 
PBOPERTIES=13 
BLOCK.  GROUP.  N1 

17  17  11 

18  18  11 

19  19  11 

C 

PAFBLOCKS  . 
ELEMENT. TYPE=36210 
PROPERTIES»14 
BLOCK.  GROUP.  N1 
20  12 

21  13 

22  12 


Y 

0 

0 

0 

0 

0 

0 


1 

-52,0 

5.0 

8 

0 

5.0 

9 

12.5 

5.0 

10 

25.0 

5,11- 

11 

37.5 

5.0 

12 

50.0 

5.0 

13 

102.0 

5.0 

14 

-52.0 

10.0 

15 

0 

10.0 

16 

12.5 

10.0 

17 

25.0 

10.0 

18 

37.5 

10.0 

19 

50.0 

10.0 

20 

102.0 

10.0 

21 

25.0 

13.0 

22 

37.5 

13.0 

23 

50.0 

13.0 

24 

102.0 

13.0 

C 

25 

76.0 

0 

26 

89.0 

0 

27 

-39.0 

5.0 

28 

-26.0 

5.0 

29 

76.0 

5.0 

30 

89.0 

5.0 

31 

-39.0 

10.0 

32 

-26.0 

10.0 

33 

76.0 

10.0 

34 

89.0 

10.0 

35 

76.0 

13,0 

36 

89.0 

13.0 

N2 

7 

7 

7 

7 

8 
8 
8 
8 
9 
9 


N2 

7 

8 
9 


N2 

7 

8 
9 


N2 

7 


20 
21 
22 

C  END  OF  PAFBLOCKS 
C 

M2SR 

reference  SPACING,  list 

1  8 

2  5 

3  § 

4  6  - 

2  I 

5  § 

9  3 


TOPOLOGY 

I, 2.  8,  9 
2,3,9.10, 

3  4,16,11 

4  5,11,12 
8,9,15,16 
9,lb,16,17 
16, li, 16, 18 

II, 12,18,19 
17,18,21,22 
18|l9,22,23 


TOPOLOGY 
5,25,12,29 
26,8,32.15 
19136, 25, 35 


TOPOLOGY 
25,26,29,30 
27; 28, 31, 32 
33  34  35,36 


TOPOLOGY 
26,  6,30,13 
7,67'14,31 

36,26,36,24 


A2 


c 
10 
11 
12 
13 
C 

C  END  OF  MESH 

MATERIAL 
^TERIAL.  NUMBER 

C 

C  END  OF  MATERIAL 
C 

PLATES .  AND .  SHELLS 


13 

5 

4,4,4,! 

4 


E 

71000 


MATERIAL . NUMBER 

11 

11 

11 


PLATE 
11 
12 

13 

14 

C 

C  END  OF  PLATES. AND. SHELLS 
RESTRAINTS 

NODE. NUMBER  PLANE 

1789  i 

C 

REPEATED . FREEDOMS 
N1  N2  PLANE 

6  0  1 

C 

C  END  OF  RESTRAINTS 
C 

SURFACE . FOR . PRESSURE 
PRESSURE.  VALUE  NODE 
-16.4  6 

C 

C  END  OF  PRESSURE 
END. OF. DATA 


NU 

0.33 


THICKNESS 

51.5 

59.5 
72 


DIRECTION 

12 

2 


DIRECTION 


PLANE 

1 


B1 


B2 


MESH 

REFERENCE  SPACING. LIST 

I  I 

1  I 

I  3 

10  13 

II  5 
C 

C  END  OP  MESH 
MATERIAL 

^TERIAL.  NUMBER  E 

J1  71 

C  END  OF  MATERIAL 
C 

PLATES . AND . SHELLS 
P^TE  MATERIAL.  NUMBER 

3-2  11 

13  11 

14  11 
C 

C  END  OP  PLATES. AND. SHELLS 
RESTRAINTS 

NODE. NUMBER  PLANE  DI 

11  1  12 

1867  1  2 

C 

C  MODELLING  OF  BOLTS 
C 

REPEATED . FREEDOMS 


E 

71000 


THICKNESS 

51 

51.5 

59.5 
72 


DIRECTION 

12 

2 


N1  N2  ■  PLANE 

6  0  1 
C 

186  501  0 

187  502  0 

8  15  0 

277  631  0 

278  632  0 

386  740  0 

387  741  0 

10  17  0 

477  831  0 

478  832  0 

C 

827  938  0 

828  939  0 

829  940  0 

830  941  0 

26  0 

920  1000  0 

921  1001  0 

C 

C  END  OF  RESTRAINTS 
C 

GAPS 

DIRECTION=2 

rn>E=i 

OFFSET=0.001 

N1  N2  NUMBER 


DIRECTION 


7  14 

173  488 

174  489 

175  490 

176  491 

177  492 

178  493 

179  494 

180  495 

181  496 

182  497 

183  498 

184  499 

185  500 

279  633 

280  634 

281  635 

282  636 

283  637 

284  638 

285  639 

286  640 

287  641 
9  16 

377  731 

378  732 

379  733 

380  734 

381  735 

382  736 

383  737 

384  738 

385  739 

479  833 

480  834 

481  835 

482  836 

483  837 

484  838 

485  839 

486  840 

487  841 
11  18 

C 


number. OF. GAP 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 


B3 


25  22  44 

820  931  45 

821  932  46 

822  933  <7 

823  934  48 

824  935  49 

825  936  50 

826  937  IS 

922  1002  52 

923  1003  53 

924  1004  54 

925  1005  55 

926  1006  56 

927  1007  57 

928  1008  58 

929  1009  59 

930  1010  60 

24  27  61° 

C 

SURFACE . FOR . PRESSDRZ 
PMSSURE.VALUE  NODE 
“16.4  6 

C 

C  END  OF  PRESSURE 
C 

END. OF. DATA 


PLANE 

1 


Cl 


APPENDIX  C:  Secondary>bending  ratio  data 

Ejqteiiroaital  secondary-bending  ratio  (SBR)  data  of  similar  specimens  from  the  United 
Kingdom  [Cl]  were  examined.  Values  of  SBR  were  obtained  from  the  finite  elonent  models 
and  con^iared  to  the  experimental  information  from  the  U.K.  specimens. 

From  [C21,  the  SBR  is  defined  as: 

STRAINbemxno 

"""  '  "  ' 

STRAINaxial 

-  (STRAIN2  -  STRAINi)  /  2 
(STRAIN2  -I-  STRAINi)  /  2 

Thus, 

STRESS2  -  STRESSi 

SBR  =  - 

STRESS2  -I-  STRESSi 

Reference  C2  also  indicates  the  position  of  the  strain  gauges  for  experimental  calculation  of 
the  SBR.  SBR  values  were  ^termined  for  the  three  finite  element  models  using  the 
infonnation  from  these  references,  where  node  positions  X  in  Figure  Cl  are  equivalent  to 
strain-gauge  locations.  Stress  values  in  adjacent  elements  to  nodes  X  were  averaged  for  the 
SBR  calculation.  Obviously  the  full  array  of  strain  gauges  could  not  be  duplicated  on  a  two- 
dimensional  model.  The  profile  of  the  specimen,  and  thus  the  profile  of  the  strain-gauge 
locations,  were  represented. 

The  values  of  SBR  for  the  three  numerical  models  remain  virtually  constant  with  change  in 
qjplied  stress  (Table  Cl)  in  contrast  to  the  experimental  data  of  reference  Cl  (Table  8)  where 
the  values  of  SBR  increase  with  increase  in  stress.  The  experimental  U.K.  values  of  SBR 
range  from  0.13, 0.22  (loading,  unloading)  at  an  applied  stress  of  17  MPa,  to  0.48,  0.42  at  an 
applied  stress  of  52  MPa,  to  0.494  at  the  maximum  applied  stress  of  104.7  MPa.  An  applied 
stress  of  104.7  MPa  is  equivalent  to  a  net-area  stress  of  350  MPa  in  the  test  sectitxi  '^f  the 
main  component.  Model  1  and  Model  2  have  rigid  connections  (repeated  frreedoms)  at  the 
location  of  the  nodes  used  in  the  SBR  calculation  which  results  in  the  negative  values  of  SBR, 
i.e.,  the  value  of  o^q)  is  unrealistically  reduced  due  to  the  excessive  restraints.  The  nodal 
stress  values  (a^d)  and  Oj(2))  increase  with  increase  in  applied  stress,  however,  they  do  so  in 
the  same  ratio,  resulting  in  ^e  same  SBR.  The  use  of  gap  elements  (as  compared  to  repeated 
freedoms)  at  the  significant  locations  of  Model  3,  improves  the  SBR  values,  i.e.,  they  are 
positive.  Model  3  predicts  the  correct  sign  of  SBR.  However,  stress  concentrations  at  the 
repeated-freedom/gap-element  interface  unduly  increase  the  value  of  thus,  the  SBR. 

The  SBR  was  determined  at  locations  either  side  of  the  strain  gauges  for  the  three  models 
(refer  to  Table  C2  and  Figure  Cl)  j  investigate  the  sensitivity  of  SBR  with  position.  The 
relatively  stiff  Model  1  has  a  snail  difference  in  SBR  values,  whereas  Models  2  and  3  display  a 
wide  range  in  values,  over  the  seven  positions,  i.e.,  the  SBR  calculations  of  the  less  rigid 
models  are  more  sensitive  to  location.  The  SBR  at  position  1  for  Model  2  could  not  be 
calculated  as  the  stress  values  in  the  adjacent  elements  were  too  extreme.  This  is  due  to 
position  1  being  the  location  of  a  repeated-freedom/gap-element  interface  for  this  model. 

The  discrepancies  between  the  experimental  and  theoretical  SBR  values  are  not  unexpected. 
The  location  of  the  strain  gauges  are  tangential  to  the  bolt-hole  edges,  hence  three- 
dimensKHial  effects  are  very  inqxvtant  Around  the  region  of  the  holes,  the  two-dimensional 
finite  dement  models  inadequacy  represent  the  strain-puge  values  of  the  actual  secotxiary- 
bending  specimen. 
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TABLE  Cl 

Secondary-bending  ratio  versus  applied  stress 
for  Model  1,  Model  2,  and  M^el  3 


Model 

No. 

Applied 

stress 

(MPa) 

Average 

stress^(u 

(MPa) 

Average 

stress^a^ 

(MPa) 

Secondary¬ 

bending 

ratio 

4.1 

9.937 

4.659 

-0.3616 

1 

16.4 

39.773 

18.647 

-0.3616 

104.7 

253.736 

118.960 

-0.3616 

4.1 

8.541 

5.482 

-0.2181 

2 

16.4 

34.213 

21.937 

-0.2186 

104.7 

218.369 

139.929 

-0.2189 

4.1 

4.383 

15.006 

0.5479 

3 

16.4 

17.660 

59.898 

0.5446 

104.7 

114.011 

380.458 

0.5389 

TABLE  C2 

Variation  in  secondary-bending  ratio  for  an  applied  stress  of  104.7  MPa 


Modd 

No. 

SBRca 

culation 

position  1 

1 

2 

3 

X 

4 

5 

6 

1 

2 

3 

-0.3424 

0.4999 

-0.3569 

-0.0413 

0.5180 

-0.3589 

-0.1419 

0.5150 

-0.3616 

-0.2189 

0.5389 

-0.3613 

-0.2543 

0.5924 

^1 

X 


representative  location  of  the  U.K.  strain  gauges 


O.  □.  ♦, 

•.  ■.  A 


location  of  pairs  of  stress  values  1  to  6  respectively 


FIGURE  Cl:  Location  of  stressx(i)  and  stressx(2)  for  the  calculation  of 

seconda^-bending  ratio 
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